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ABSTRACT 

We re-examine the photometric and spectroscopic evidence available for the star OGLE-TR-56, 
recently discovered to harbor a giant planet presenting transits and orbiting with a period of 1.21 
days. We investigate the possibility that the observational signatures reported might be the result 
of "blending" with the light of an eclipsing binary along the same line of sight. Using techniques 
developed earlier we perform fits to the light curve under a variety of blend configurations, subject 
to all observational constraints, and we infer further properties of these possible blends. We then 
carry out realistic end-to-end simulations based on those properties in order to quantify the spectral 
line asymmetries and radial velocity variations expected from such scenarios. We confront these 
calculations with the observations. The results from these models are clearly inconsistent with the 
measured radial velocity and bisector span variations, ruling out blends and confirming the planetary 
nature of the companion. The example of OGLE-TR-56 serves to illustrate the sort of tests that 
can and should be performed on transiting planet candidates to eliminate the possibility of "false 
positives" , and in particular, line-of-sight contamination from eclipsing binaries. 
Subject headings: binaries: eclipsing — line: profiles — planetary systems — stars: evolution — stars: 
individual (OGLE-TR-56) — techniques: radial velocities 



1. INTRODUCTION 

The search for extrasolar giant planets has entered 
an exciting new era, one in which these objects can 
be discovered indirectly by the tiny photometric signa- 
tures they produce as they transit across the disk of 
their parent star (periodic ~1% drops in brightness). 
Numerous photometric search programs are underwa y 
or under development worldwide (see iHornei |2003). 
Six stars with t ransiting planets have been found so 
far: HP 209458 iHenrv et alJ 120001: jCharbonneau et"aD 
l2000h OGLE-TR-56 llKonacki et aljboOSal). O GLE-TR- 
113 llBouchv et all l2004UKonacki et alJl200l. OGLE- 
TR-1 32 ijBouchv et alJ 120041). TrES^l llAlonso et all 
2004), and OOLE-TR-1 1 1 ijPont et aTl 120041) 4 The 
OGLE-TR planets were originally detected as low- 
amplitude transit events by the most successful of the 
photometric surveys, OGLE-III (Optical Gravita tional 
Lensing Experiment: lUdalski et aLll2002blcL 120031) . and 
later confirmed by radial-velocity follow-up. The planet 
around HD 209458 was originally discovered by its radial 
velocity signature — the small reflex motion of the star in 
response to the pull from the planet — and only later was 
it found to undergo transits. It is the Doppler technique, 
in fact, that has been the most successful in discovering 
extrasolar plane ts, having produ ced more than 100 giant 
planets to date l)Schneiderll2004|) . 

Both the Doppler and the transit search methods suf- 
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for from the problem of astrophysical "false positives", 
which refers to other phenomena that can produce the 
same observational signature but do not involve plan- 
ets. In the case of the Doppler searches the concern is 
that similar velocity variations might result from chromo- 
sphcric activity on the star (e.g., spots), sub-surface con- 
vection (granulation), other line asymmetries from un- 
detected stellar companions, or even stellar oscillation s 
(see, e.g.. lGravlll"997t iGrav fc HatzesllT997t lGravlll998j) . 
For the transit searches the main contaminants are (1) 
eclipsing binaries with grazing geometry, (2) transits of 
a small star in front of a large star (such as an M dwarf 
passing in front of an F star, or a solar-type star in front 
of a giant), and (3) eclipsing binary systems with deep 
eclipses, blended with the light of another star along the 
same line of sight that dilutes the minima. In the latter 
case (referred to as a "blend"), the intruding star may or 
may not be physically associated with the binary. It is 
these blend scenarios that are the subject of this paper. 

Cases that only involve an eclipsing binary with no con- 
taminants (scenarios [1] and [2] above) are relatively easy 
to dismiss with radial velocity measurements since the 
expected amplitudes are typically several tens of kms -1 , 
as opposed to a few hundred ms _1 for the case of a 
Jupiter-mass planet around a solar-type star. Blend 
configurations, on the other hand, can be much more 
difficult to rule out as the main star may show little 
or no velocity variation at all, and may therefore ap- 
pear to have a very low-mass companion in orbit when 
in reality it does not. This is particularly the case 
for very faint candidates (V = 14-19) such as those 
being reported from surveys suc h as OGLE-III (e.g., 
l Udalski et ai]l2002aD. EXPLORE <|Mallen-Ornelas et alJ 
I2003D . MACHO iforake Cookll2004l) . and others. For 
these objects the precision of the velocities is much lower 
than for the typically brighter stars (V < 11) that have 
been targeted in Doppler searches, and the potential for 
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confusion is greater because they also tend to be in very 
crowded fields. It is therefore important to examine these 
cases very carefully. 

In a recent paper (jTorres et al.ll2004bl hereafter Pa- 
per I) we described a procedure to model the light curves 
of transit candidates as the combination of three stars 
forming the blend, and we applied it to the case of 
OGLE-TR-33, a candidate that turned out to be a false 
positive. We showed how such modeling tools allow one 
to predict other observable properties of the blend and 
can become a useful guide to look for additional signs of 
a contaminating eclipsing binary. One such sign that can 
easily be overlooked is asymmetries in the spectral lines, 
which can lead to spurious velocity measurements. In 
the present paper we develop procedures to simulate line 
asymmetries and their effects on the measured velocities 
of a transit candidate resulting from typical blend scenar- 
ios. We illustrate these techniques by carrying out a crit- 
ical examination of the case of OGLE-TR-56 (V — 16.6), 
the first extrasolar transiting planet to emerge from the 
lists of very faint candidates, and the one with the short- 
est orbital period so far (1.21 days). We perform blend 
simulations subject to all observational constraints in an 
attempt to explain the system as a false positive, but 
despite all our efforts we are unable to find a plausible 
configuration involving a contaminating eclipsing binary 
that does not violate one or another of those constraints 
by a large margin. This lends further support to the 
planetary nature of the companion of OGLE-TR-56. 

2. OGLE-TR-56: THE SPECTROSCOPIC EVIDENCE 

The original detection of peri odic transit-like event s 
in OGLE-TR-56 was reported bv lUdalski etaD ((2002b). 
The star was part of the second batch of candidates 
produced by the OGLE-III photometric survey, which 
targeted three fields toward the Galactic center. High- 
resolution spectroscopic follow-up observations were con- 
ducted by our team in July 2002, and small but sig- 
nificant radial velocit y variations were announced by 
IKonacki et al.l l|2003a|) on the basis of measurements on 
three nights near the quadratures, implying the presence 
of a planetary-mass object in orbit. As the first such 
claim among all the transit candidates reported until 
then — and by far the faintest candidate ever considered 
to have a chance of harboring a planet — those results 
were carefully scrutinized by the community and were 
received quite understandably with a measure of skepti- 
cism. The initial concerns on the observational side re- 
garding this particular detection focussed on two issues: 
(a) The reality (or statistical significance) of the radial 
velocity variations, especially since the claim was based 
on only 3 averaged measurements; and (b) the interpre- 
tation of the velocity variations, given the crowding in 
the field and the increased probability of a blend. 5 

Such concerns are likely to be repeated for other faint 

5 Two additional concerns have to do with (c) whether a planet in 
such a tight orbit such as OGLE-TR- 56 can exist at all (which the- 
ory appears not to rule out; sec Trilling ct al. 1998; Baraffc ct al. 
2004) , and (d) why Doppler searches among brighter stars have not 
found additional examples in very short-period orbits (< 2.5 days) 
despite being sensitive to them. The latter may perhaps be the re- 
sult of a much lower frequency of occurrence of these "very hot 
Jupiters" , coupled with very strong selectio n toward short pe- 
riods in transit surveys such as OGL E (see IBouchv' et ail 2004: 
IGaudi. Seager. fc Mallen-Orneladl200l) . 



candidates that might be confirmed in the future, and 
it is therefore necessary to consider them very seriously. 
The first concern is related mostly to errors of measure- 
ment (accidental or systematic), and whether the re- 
sults might be due, for example, to limitations in the 
reduction and/or analysis techniques or to the stability 
of the spectrograph employed (in thi s case, HIRES on 
the Keck I telescope; iVoet et alj[l99l. For OGLE-TR - 
56 these issues were addressed by IKonacki et al.l (|2p03b) , 
who showed that systematic errors in the velocities are 
no more than about 100 ms" 1 and accidental errors are 
somewhat smaller, making the reported peak-to-peak ve- 
locity change (initially estimated to be ~ 340 ms" 1 ) sta- 
tistically significant considering that the phasing in the 
orbit is known and fixed from the photometry. Eight ad- 
ditional measurements w ith the same instrum ent were re- 
ported more recently bv lTorres et alJ l|2004aj) . and show 
even more clearly that the velocity is changing as ex- 
pected from the phasing in the orbit. The revised ampli- 
tude folding together all the measurements is larger than 
originally estimated, but it is now much better charac- 
terized (peak-to-peak variation near 600 ms -1 ). More 
importantly, independent confirmation of velocity vari- 
ations has been repor ted by another grou p using differ- 
ent instrumentation. iMavor et all l)2003|) obtained five 
measurements of OGLE-TR-56 during the commission- 
ing of the new HARPS instrument on the ESO 3.6-m 
telescope in June of 2003, achieving a nominal precision 
of about 30 ms -1 . Thei r preliminary analys is gives much 
the same amplitude as lTorres et al.l (|2004a|K 

With the reality of the velocity variations now beyond 
doubt, the issue of their interpretation becomes the main 
concern for OGLE-TR-56, and is the motivation for the 
remainder of this paper. Can the observations be ex- 
plained by a blend scenario, with an eclipsing binary 
along the same line of sight contaminating the light of 
the main star? If so, the velocity variations measured 
would be caused by blending with the spectral lines from 
the primary of the eclipsing binary moving back and forth 
relative to the main star, with a period of 1.21 days. Are 
these extra lines visible in the spectra? Strong and vari- 
able line asymmetries are also to be expect ed from such 
a scen ario. Are the measured asymmetries ({Torres et alJ 
l2004a|) consistent with a blend model? Given the impor- 
tance of OGLE-TR-56 as the first example of a seemingly 
new class of planetary objects under extreme conditions 
of proximity to the parent star, the answers to the ques- 
tions above must be quantitative rather than qualitative. 
We begin in the next section by modeling the light curve 
in detail following the procedures we introduced in Pa- 
per I for OGLE-TR-33. 

3. LIGHT CURVE FITS TO OGLE-TR-56 UNDER A BLEND 
SCENARIO 

The shape of a light curve showing transit-like events 
contains valuable information that can be used in some 
cas es to recogniz e fal se positive scenarios. Fo r exam- 
ple. iDrakd l|200l and lSirko fc Paczvhskil lp¥)3|) showed 
that the variations outside of eclipse, as measured from 
a Fourier analysis, can often indicate that the compan- 
ion is stellar as opposed to substellar, from effects due to 
tidal distortions, reflection, or gravity brightening. An 
isolated eclipsing binary (e.g., one with grazing geome- 
try, or consisting of a small star passing in front of a 
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much larger star) can be ruled out as the explanation for 
the light curve of OGLE-TR-56 because of the very small 
velocity amplitude measured. However, an eclipsing bi- 
nary blended with a another star remains a possibility, 
in principle, and some of its properties could still show 
through in the light curve, even though they would be di- 
luted. But as reported by the authors above, the Fourier 
decomposition test for OGLE-TR-56 comes out negative, 
indicating no significant curvature outside of eclipse. 

A different kind o f test was proposed by 
iSeaeer fc Mallen-Ornelasl l)2003f) based on the mor- 
phology of the transits themselves. They showed that 
the mean density of the primary star, as well as the 
radius of the companion, can be estimated by careful 
measurement of the depth of the transit (relative change 
in flux, AF), its total duration (tx), and the duration 
of the flat portion (tp), along with the knowledge of 
the orbital period. Any inconsistency in the mean 
density with values for normal (main-sequence) stars, or 
with additional information available for the primary, 
would be an indication of a blend. Once again this 
test comes out negative for OGLE-TR-56. We measure 
AF = 0.016, t T = 0.060, and t F = 0.036 (the latter two 
in phase units), which w hen combined with t he orbital 
period of 1.2119189 days ijTorres et alJl2004aj) lead to a 
primary star very similar to the Sun in mean density, 
as well as in both mass and radius when a typical 
mass-radius relation for main sequence stars is adopted. 
Indeed the temperature estimated for OGLE-TR-56 
from our high -resolution spectra is solar-like (~5900 K; 
!Sasselovll2003l) . 

More detailed information on the properties of a pos- 
sible blend can be obtained from a direct fit to the light 
curve as described in Paper I. Briefly, the photometry is 
mod eled using the eclipsing binary computer code EBOP 
(see [Nelson fc David 11973 lEtzei Il98ll iPopper fc Etzel 
1981), fully accounting for effects such as stellar oblate- 
ness, reflection, limb darkening, and gravity brightening. 
The properties of the three stars composing the blend 
(which are assumed to be on the main sequence) are pa- 
rameterized in terms of their mass and are constrained 
to fit suitable model isochrones. These isochrones may 
be different for the binary and for the main star if they 
do not form a physical triple system. In that case, differ- 
ential extinction must also be accounted for since they 
can be at different distances from the observer. All other 
relevant properties of the stars such as their radii and lu- 
minosities are taken also from the isochrones (see Paper I 
for details). The important advantage of this procedure 
is that it allows testable predictions for some of the other 
properties of the eclipsing binary that are potentially ob- 
servable. These properties include the brightness of the 
primary in the eclipsing binary relative to the main star, 
its velocity semi-amplitude, and its projected rotational 
velocity (assuming its spin is synchronized with the or- 
bital motion). This was illustrated in Paper I with the 
example of the blend case OGLE-TR-33. The same pro- 
cedure is applied here to OGLE-TR-56. 

As in Paper I we m ake use of the model isochrones 
bv lGirardi et alJ (J2000) to describe the properties of the 
stars in the blend 6 . The photometric data for OGLE-TR- 

6 For low-mass stars these models become increasingly unreal- 
istic partly because of shortcomings in the treatment of convec- 



56 are described bv lTorres et alJ l|2004aj) . and we adopt 
the measurement errors reported by the OGLE team 7 . 
Those observations were made in the I band. We be- 
gin by assuming the eclipsing binary and the third star 
are physically associated, and we therefore use the same 
isoch rone for both with an adopted age of 3 Gyr ijSasselovl 
12003ft and solar metallicity (Z = 0.019). We refer to this 
model as MODEL 1. We follow the notation of Paper I and 
refer to the primary and secondary in the eclipsing binary 
as star 1 and star 2, respectively, while the main star that 
dilutes the eclipses of the binary is referred to as star 3. 
We assign to star 3 a mass on the isochrone such that 
the corresponding temperature agrees with the spectro- 
scopic determination for OGLE-TR-56 (M 3 = 1.05 M , 
T ff' = 5900 K). The inclination angle is assumed to be 
90° for simplicity. 
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Fig. 1.— (a) Blend model fit to the light curve of OGLE-TR-56 
(M0DEL1). The main star (star 3) and the contaminating eclipsing 
binary are assumed to be physically associated (age = 3 Gyr). The 
masses of the stars in the eclipsing binary (star 1 and star 2) are 
Mi = 0.89 M and M 2 = 0.13 M , and star 3 contributes 71% 
of the light in the I band. The predicted light ratio in the visual 
band is L1/L3 = 0.35. (b) Enlargement of the primary eclipse. 

tion, in the opacities, in the equation of state, and the use of gray 
boundary conditions, all of which become very important in the 
lower main sequence. In particular, theoretical radii from these 
and all current stellar evolution models are known to be underesti- 
mated by 10% or more compared to measurements of well-s tudied 
low-m ass eclipsing binaries (see, e.g., Torres & Ribas 2002; Ribas 
120031) . While ot her models for cool stars are somewhat more so- 
phist icated (e.g.. I^Tess. Forestini fcTTo ueados 1997; Baraf fe~et alJ 
1998), our reason for choosing this particular set of isochrones is 
that they span the largest range of masses for our application, be- 
ginning at 0.15 Mq. Despite this, slight extrapolations to even 
lower masses were necessary in some cases. Given the importance 
of the model radii in our simulations for determining key properties 
of the eclipses such as their depth and duration, we have attempted 
to bring them into closer agreement with the radii of real stars by 
applying ad-hoc corrections based on a careful comparison with 
accurately measured radii for late-type stars. The correction fac- 
tors are mass-dependent, and average 1.1 over the mass range of 
interest for our work. 

7 See http: / /bulge .princeton. edu/~ ogle/ogle3/transits/ 
ogle56.html, and Kruszcwski & Scmcniuk 12003). 
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Our blend model adjusted to the light curve is shown 
in Figure H The resulting properties of the stars in the 
eclipsing binary are M\ — 0.89 M Q (approximately spec- 
tral type G8 V) and M 2 = 0.13 M Q (M4-5 V), and the 
inferred distance to the system is ~1.8 kpc. Very shallow 
secondary eclipses only 0.002 mag deep are also predicted 
by this model, but they are undetectable given the pre- 
cision of the observations (typical errors of ^0.006 mag 
for an individual measurement). We note that the fit to 
the primary eclipse with this blend model is essentially 
indistinguis hable from a fit that models a true planetary 
transit (see iTorres et al.l r2004al Figure 4). Thus, the 
shape of the light curve alone is insufficient to discrim- 
inate between the two. The fitting procedure provides 
all the physical properties of star 1 from the isochrone, 
so that it is possible to estimate the semi-amplitude of 
its radial- velocity orbit (K\ — 26 kins" 1 ) and also its 
projected rotational velocity (v% sirxi — 36 kms -1 ), un- 
der the plausible assumption that it is synchronized with 
the orbital motion (given the short orbital period). The 
brightness of star 1 relative to the main star is predicted 
to be Li/L 3 = 0.35 in the visual band 8 . The latter 
quantity conflicts with the spectroscopic observations for 
OGLE-TR-56, which place an upper limit on the pres- 
ence of lines from another star in the spectrum of ab out 
3% of the light of the main star ( Konacki et al. 2003b) . A 
star as bright as 35% of the light of star 3 would be fairly 
obvious in the spectra, despite the rotational broaden- 
ing. We conclude that this particular blend scenario is 
inconsistent with the data 9 . Additional difficulties with 
MODEL 1 arise from the magnitude of the line asymmetries 
that are expected. We defer the discussion of these issues 
to gl 

Since the main problem with MODEL 1 seems to be the 
excessive brightness of star 1, we explored ways of mak- 
ing star 1 fainter while still providing a reasonably good 
fit to the light curve. One such way is to remove the 
requirement that the eclipsing binary and the main star 
be in a physical triple, and to allow the binary to be in 
the background. In this case the age of the isochrone for 
the binary can in principle be different. For lack of a 
better estimate we retain the 3 Gyr isochrone for both, 
but the binary is no longer constrained to be at the same 
distance as star 3. We performed extensive modeling of 
the photometry for a range of distances between the bi- 
nary and star 3 parameterized, for convenience, in terms 
of the difference in distance modulus AM (see Paper I). 
We assumed an edge-on orientation, as before. The re- 
sults are summarized in Figure^ where we show how the 
parameters of star 1 and star 2 change as the binary is 
placed farther behind star 3. As expected the observed 
light ratio L1/L3 decreases with distance (Figure^), 
and reaches the upper limit determined spcctroscopically 

8 The symbol L is used throughout this paper to denote bright- 
ness ("light"). 

9 Tests show that adjusting the age of the isochrone does not 
qualitatively change the conclusion. With increasing age L1/L3 
decreases somewhat as star 3 evolves upward in luminosity, but not 
enough before a star with a mass appropriate for the temperature 
of OGLE-TR-56 leaves the main sequence altogether. The latter 
occurs for ages greater than about 7 Gyr, at which L1/L3 has 
only decreased to approximately 0.20. Further tests with lower 
inclination angles for the binary show that the L1/L3 increases 
as the system departs from an edge-on orientation, making the 
disagreement with the observations even worse. 



(0.03) at AM w 2.0 (binary approximately 2.3 kpc be- 
hind star 3; Figure Et). We refer to this solution as 
M0DEL2, and indicate it with squares in the figure. The 
mass of star 1 that provides the best fit has increased at 
this distance by about 12% compared to M0DEL1, while 
that of star 2 is more than a factor of 3 larger than be- 
fore (Figure |3i) . The reason for these changes can be 
understood as follows. As the binary is pushed back and 
its contribution to the total light decreases, the eclipse 
becomes more diluted (shallower), and so in order to 
maintain the observed depth the relative size of the sec- 
ondary in the binary (star 2) must increase relative to 
the primary (star 1), which implies increasing the mass 
of star 2. This, in turn, will affect the total duration 
of the eclipse for a fixed star 1, both from the change 
in the relative size of star 2 and star 1, and from the 
change in semimajor axis for a fixed orbital period (see 
iSeaeer fc Mallen-Ornela s 2003). In order to conform to 
the observed transit duration, the mass and size of star 1 
then need to increase slightly to compensate. Accord- 
ingly, the predicted v sin i of star 1 shows only a small 
dependence with AM, while K\ varies significantly as 
the mass of the companion increases. The quality of the 
fits degrades as the binary is placed farther in the back- 
ground (Figure Hf), but remains more or less acceptable 
at AM = 2.0. That fit is shown later in Figure^ along 
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Fig. 2. — Results from a grid of light curve fits to the photome- 
try of OGLE-TR-56, with the eclipsing binary placed at increasing 
distances behind the main star (star 3), as parameterized by the dif- 
ference in the distance modulus AAf. The fit marked with squares 
is referred to in the text and in Table 1 as M0DEL2. The leftmost 
points correspond to M0DEL1. (a) Fitted distance in parsecs to the 
binary and to star 3; (b) Light ratio L1/L3 in the visual band. 
The horizontal dotted line represents the upper limit estimated 
from the high-resolution spectra (~3%); (c) Contribution of star 3 
to the total light in the I band ("third light" in eclipsing binary 
terminology); (d) Best fit masses for the primary and secondary in 
the eclipsing binary; (e) Predicted radial velocity semi-amplitude 
and projected rotational velocity for the primary in the eclipsing 
binary (star 1); (f) Reduced \ 2 for the solution (that the values 
are lower than unity suggests the photometric errors are slightly 
overestimated) . 
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with the light curve from MODEL 1 for comparison. The 
consequences of MDDEL2 for the line asymmetries and ex- 
pected velocity variations are explored in the next sec- 
tion. 
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Fig. 3 . — Isochrones from the Yo nsei-Yale series of evolutionary 
models ( Yi. Kim, & Dcmarquc 2003) for a range of ages labeled in 
Gyr. The 3 Gyr isochrone for OGLE-TR-56 is represented with a 
heavy dashed line. As an example, the asterisks show two possible 
locations of a main-sequence star at a fixed age with the same ef- 
fective temperature, but with brightness differing by about 2 mag. 
The effect of choosing the brighter location for star 3 in a blend 
model would be to decrease L1/L3. The open circle on the 3-Gyr 
isochrone corresponds to the position of OGLE-TR-56 on the main 
sequence, at the measured temperature of 5900 K (vertical dotted 
line). At these low temperatures the flexibility of choosing the 
brighter of two vertically aligned locations for star 3 with signifi- 
cantly different brightness is effectively eliminated because of the 
change in the shape of the isochrones. 

If we wish to keep the three stars at the same distance, 
as in a triple system, an alternate way of decreasing the 
relative brightness of star 1 to comply with the spec- 
troscopic constraint is to make star 3 itself intrinsically 
brighter, e.g., by evolving it from the zero-age main se- 
quence. There are good reasons for assuming physical 
association that have to do with the apparent symmetry 
of the velocity variations (see <0} . The brightness change 
required here is roughly a factor of 10 (about 2.5 mag- 
nitudes) in order to reduce the brightness ratio L1/L3 
from 0.35 (MODEL 1) to the observational upper limit of 
0.03 ijKonacki et al. 2003b). A similar situation was de- 
scribed in Paper I for the case of OGLE-TR-33, where it 
was shown that by fine-tuning the age an isochrone could 
be found in which a slightly more evolved star 3 could be 
about 2 magnitudes brighter at the same effective tem- 
perature (constrained by the spectroscopy). The location 
of star 3 in the H-R diagram would be beyond the turnoff 
but still on the main sequence, n ear the point of hy - 
drogen exhaustion (see Figure 5 of lTorres et alJ l2004b). 
OGLE-TR-56 is cooler than OGLE-TR-33, however, and 
a glance at the morphology of the isochrones reveals that 
a similar argument does not work for a temperature of 



5900 K, at an y age. This is i l lustra ted in Figure^ where 
instead of the iGirardi et all 1120001) m odels we have used 
those of lYi. Kim, fc Demaraud l)2003|) for their finer res- 
olution near the turnoff. It is seen, for example, that 
the "blue hook" feature gradually disappears for older 
ages, reducing the difference in brightness between its 
red end (coolest and most luminous point on the main 
sequence beyond the turnoff) and a location below it with 
the same temperature. The change in shape as the age 
increases coincides with the disappearance of convective 
cores at masses slightly larger than the Sun. Unless one 
is willing to accept that OGLE-TR-56 is in the rapid 
phase of evolution often referred to as the Hertzprung 
gap (horizontal portions of the isochrones where stars 
are evolving mostly toward cooler temperatures without 
changing their brightness significantly), the brightness of 
star 3 cannot be much larger than indicated by the open 
circle, at its location on the main sequence. 
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Fig. 4.— Two blend models for OGLE-TR-56 (thick lines) that 
satisfy the spectroscopic constraint on the brightness of star 1 rel- 
ative to star 3 in the V band (< 3%). The fit for M0DEL1 from 
Figure ^> is superimposed in each panel for reference (thin line), 
(a) M0DEL2, in which the binary is in the background of star 3 
(~2.3 kpc behind, in this case; see Figurel2l, (b) M0DEL3, in which 
the binary and star 3 are at the same distance (as in a physical 
triple), but star f is constrained to have a brightness of 3% relative 
to the main star. 

For completeness we explored also a blend scenario in 
which the three stars form a physical triple system, but 
star 1 is constrained to have a relative brightness in the V 
band consistent with the upper limit from spectroscopy 
(~3%). Only the properties of star 2 were adjusted, and 
we refer to this as M0DEL3. The inclination angle was 
again fixed at 90°, and the adopted age is the same as be- 
fore (3 Gyr), although this does not qualitatively change 
the results. This solution is noticeably worse than the 
previous two (reduced \ 2 about 8% larger), but at least 
it still reproduces the total duration and depth of the 
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TABLE 1 



PcLr&niotcr 


MODEL1 


MODEL2 


MODEL3 


Mi (Mq) 


0.89 


1.00 


0.60 


M 2 (Mq) 


0.13 


0.42 


0.22 


L1/L3 (V band) 


0.35 


0.03 


0.03 


L 3 /L to t (I band) (%) a 


71 


94 


94 


K x (kms" 1 ) 


26 


66 


50 


vi s'mi (kms -1 ) 


36 


40 


26 


Distance to star 3 (kpc) 


1.78 


1.54 


1.54 


Distance to binary (kpc) . . . 


1.78 


3.86 


1.54 


AM (mag) b 


0.00 


2.00 


0.00 


a v (mag kpc _1 ) c 


0.56 


0.66 


0.66 


Reduced \ 2 of the n t d 


0.78188 


0.78609 


0.84734 



sians with appropriate parameters to simulate the mea- 
sured correlation functions. From numerical experiments 
they were able to infer the key parameters of the system, 
and to fully explain the observations as the result of a 
possible brown dwarf orbiting the visual companion in- 
stead of a planet around the primary star. Their results 
for the physical and orbital properties of th e visual com- 
panion were subsequently corroborated by iZucker et alJ 

In this paper the situation is somewhat different in 
that our blend models fitted to the light curve of OGLE- 
TR-56 constrain many of the key properties of the stars. 
Our goal is then to see if these inferred properties lead 
to a consistent picture of velocity and asymmetry vari- 
ations matching what is observed. In addition, we wish 
to be as realistic as possible in approximating the obser- 
vations. Therefore, instead of adopting Gaussian profiles 
for the correlation functions, we simulate composite spec- 
tra guided by the models in JJ3 by combining together 
calculated spectra for star 1 and star 3 that conform 
to the stellar and orbital characteristics implied by each 
blend scenario. We then measure radial velocities by 
cross-correlation and estimate spectral line asymmetries 
in these artificial spectra in exactly the same way as was 
done for the original observations of OGLE-TR-56. 

For star 3 the calculated spectrum we used is the same 
one adopted bv lTorres et al.1 l)2004al) . which was derived 
from a fit to the observed spectrum of OGLE-TR-56, 
and was used in that work as the template for the cross- 
correlations. This calculated spectrum has an effective 
temperature of 6000 K, surface gravity logg = 4.5, and 
projected rotational velocity usini = 3 kms -1 . For a 
given blend model star 1 is assigned the temperature and 
surface gravity inferred from the isochrone for its mass, 
and a rotational velocity given by the predicted value 
of v\ sini. Star 1 was then scaled down by the predicted 
light ratio L1/L3, and Doppler-shifted before adding it to 
the spectrum of the main star. The shift between star 1 
and star 3 is controlled by the velocity amplitude of star 1 
in its orbit (-Ki) and the difference A7 between its center- 
of-mass velocity and the (presumably constant) velocity 
of star 3. Since the period and epoch of transit are well 
known, the Doppler shift can be computed for any instant 
of time. A7 is a free parameter in principle, although for 
scenarios such as MODEL 1 and M0DEL3 it is expected to 
be small compared to K\ given that the configuration is 
assumed to be a hierarchical triple system. 

In an effort to approximate the real observations of 
OGLE-TR-56 as closely as possible, the artificial compos- 
ite spectra were divided into sections corresponding pre- 
cisely to e ach of the echelle orde rs of the original observa- 
tions (see lKonacki et a~ 2003b), and rebinned to exactly 
the same number of pixels as in the real data prior to fur- 
ther analysis. The calculated spectra were initially gen- 
erated at much higher spectral resolution than the real 
observations (by a factor of ~10), and then broadened 
with the average instrumental profile of the Keck spec- 
trograph, which was assumed to be Gaussian in shape 
and with FWHM = 4.6 kms -1 . 

Radial velocities from these artificial composite spec- 
tra were measured by cross-correlation using the IRAF 10 

10 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities 



a Relative contribution of star 3 to the total light of the system. 
b Difference in the distance modulus between the binary and 
star 3. 

c CoefBcient of differential extinction in the V band, adjusted to 
reproduce the observed system magnitudes in the V and I bands 
(see Paper I). 

d The values significantly lower than unity imply that the mea- 
surement errors are slightly overestimated, by roughly 10%. 



eclipse tolerably well (see Figure 0Jd). We summarize the 
properties of this and the other two models in Table 1. 

As in the case of the first scenario, M0DEL2 and 
M0DEL3 predict shallow secondary eclipses with depths of 
^0.002 mag and ^0.006 mag, respectively, which would 
be difficult to see. 

4. SIMULATING LINE ASYMMETRIES AND RADIAL 
VELOCITY VARIATIONS 

In blend scenarios such as those modeled in the previ- 
ous section the presence of lines from another star in the 
spectrum of OGLE-TR-56 will affect the measurement of 
the radial velocities to some degree. As mentioned above 
this is due to asymmetries introduced in the profiles of 
the main star. The extent to which these asymmetries 
will distort the radial velocities depends on the brightness 
of the secondary lines and their relative displacement. 
Intrinsic broadening may also be a factor. Quantitative 
measures of the asymmetries for planet candidates show- 
ing velocity variations are now often reported in the lit- 
erature, and in many cases they have turned out to be 
small. We believe, however, that this alone is not suffi- 
cient evidence against a blend. In this section we explore 
the magnitude of these effects by means of simulations to 
see whether they can explain the measure ments reported 
for OGLE-TR-56 bv lTorres et all ll2004al) . 

A s imilar philosophy was adopted by iSantos et alJ 
( 2002) to investigate small radial velocity variations and 
line asymmetries found in the case of the star HD 41004. 
This is a close visual binary initially thought to have a 
planetary companion around the primary star based on 
high-precision Doppler monitoring performed by those 
authors. The visual binary with an angular separation of 
0'.'5 was unresolved at the spectrograph, effectively mak- 
ing it a blend. Unlike the case of OGLE-TR-56, however, 
there is no transit signature in HD 41004 and therefore 
there is less infor mation about the sta rs and the geome- 
try of the system. ISantos et al.l l)2002t) approximated the 
cross-correlation functions of the two stars in HD 41004 
with Gaussians, and sought to reproduce the measured 
velocity and bisector variations by combining two Gaus- 
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task XCSAO ijKurtz fc Minklll99Sj) . combining the re- 
sults for all orders. Th e template used is th e same as for 
the real observations ijTorres et alJl2004a|) . and it hap- 
pens to be the same spectrum used to model star 3. Line 
asymmetr ies were qua ntified by computing the line bisec- 
tors (e.g.. IGravlll992|) directly from the co-added corre- 
lation function, which is representative of the average 
spectral line. We then calculated the velocity difference 
in the bisector at two different correlation levels, which 
we refer to as the "bisector sp an". 11 This i s anal ogous 
to the defini t ion em ployed by iQueloz et alJ 1)20011) and 
ISantos et all (|2002j) . 

We first simulated MODEL 1, sampling all phases of the 
binary orbit and assuming A7 = kms -1 . The results 
are shown in Figure |SJ and are compared with the ra- 
dial velocity obser vations and b i sector spans reported for 
OGLE-TR-56 bv iTorres et all l|2004a(l . The top panel 
indicates an excellent agreement with the measured ve- 
locities, which is however only a coincidence. The lower 
panel shows that the asymmetries expected from this sce- 
nario are much larger than observed. The predicted bi- 
sector spans can reach velocities exceeding ±500 ms _1 
near the quadratures, while the observations show no sig- 
nificant variation at all within the errors. Thus M0DEL1 
fails to reproduce the observations, as was already an- 

for Research in Astronomy, Inc., under contract with the National 
Science Foundation. 

11 The particular correlation levels used here, 0.2 and 0.4, were 
chosen from the properties of the correlation functions for the ob- 
servations of OGLE-TR-56. Levels below 0.2 are affected by noise, 
and levels above 0.4 are not measurable for some of the weaker 
exposures in which the peak of the correlation function does not 
reach that value. Our conclusions below are qualitatively unaf- 
fected by chosing different levels between these two values, so long 
as the simulations use the same levels. 
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ticipated in fJ3 from the excessive brightness of star 1 
relative to star 3 (Li/L 3 = 0.35). 

Next we explored the sensitivity of the radial velocity 
and bisector span variations to the main parameters of 
the blend, in an attempt to understand their effect. We 
carried out extensive simulations over a broad range of 
values of L1/L3, K±, A7, and Disini, varying one pa- 
rameter at a time and leaving the others fixed at their 
nominal values for MODEL 1. A sampling of these tests is 
described below. 

The effect of imposing different ratios L1/L3 for the 
brightness of the primary of the eclipsing binary relative 
to the main star in the visual band is shown in Figure El 
As expected the amplitude of the variations scales with 
the brightness of star 1, although not quite linearly. 

Changing the semi-amplitude K\ results in relatively 
minor changes in the amplitudes, but significant changes 
in the shape of the simulated curves both for the bisector 
span and the velocities, as seen in Figure [7| For small 
values of K\ the lines of star 1 remain within the wings 
of those of the main star at all phases, and the curves are 
nearly sinusoidal. As K\ increases, the lines of star 1 shift 
out of the wings of star 3 at the quadratures, and the ef- 
fect is smaller at those phases producing a flattening out 
of the curves. They beco me increasing l y non -sinusoidal 
with larger values of K\ . ISantos et alJ l|2002fl described 
the same effect in their simulations. 

Both the amplitude and the shape of the predicted 
velocity variations and bisector span variations depend 
very sensitively on the rotational broadening of star 1. 
This is illustrated in Figure [5J Sharper lines for star 1 
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Fig. 7. — Effect of varying the semi-amplitude K\ (in kms x ) 
upon the expected radial velocity variations (left panels) and bi- 
sector spans (right panels) from a blend scenario. K\ = 26 kms -1 
corresponds to MQDEL1, and the remaining parameters of the simu- 
lations are held at their values in that model (see Table 1). 
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Fig. 9.— Effect of A7 (in kms -1 ) on the expected radial veloc- 
ity variations (left panels) and bisector spans (right panels) from 
a blend scenario. M0DEL1 has A7 = kms -1 . The remaining pa- 
rameters of the simulations are held at their values in that model 
(see Table 1). 
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Fig. 8. — Effect of nisini (in kms -1 ) on the expected radial 
velocity variations (left panels) and bisector spans (right panels) 
from a blend scenario. Note the change in the vertical scale for the 
bottom two panels on each side (see text). MDDEL1 has t;isini = 
36 kms -1 . The remaining parameters of the simulations are held 
at their values in that model (see Table 1). 



give larger variations except near the quadratures (for 
K\ = 26 kms -1 , as in M0DEL1), as might also be ex- 
pected. This leads to cusps in the curves that are pro- 
gressively closer to phase 0.0 and 0.5 as the rotational 
velocity becomes smaller, qualitatively resembling the 
effect of a larger K\ value. As vishu increases, the 
effect is attenuated due to the broader lines and the 
curves approach a more sinusoidal shape. In particular, 
a large drop in amplitude is observed between v\ sin i of 
20 kms -1 and 30 kms -1 , coincident with the point at 
which the rotation numerically exceeds the value of K\. 

Finally, a change in A7 produces another qualitative 
distinction, this time an asymmetric shape to the bisec- 
tor span and velocity variations that depends upon the 
sign of A7. This is shown in Figure ED As a result, 
blend models that do not assume physical association 
between the eclipsing binary and the main star (such as 
M0DEL2) would be expected to show asymmetrical shapes 
to the bisector and velocity variations, since A7 can be 
arbitrarily large in those cases. The measured velocity 
variations for OGLE-TR-56 appear to be more or less 
sinusoidal (c.f. Figure EJ), and therefore values of A7 
significantly different from zero (of order 10 kms -1 or 
greater, given the number of observations and their un- 
certainties) would seem to be inconsistent with the ob- 
servations. Similarly, values of K\ much larger than that 
used in MODEL 1 would lead to a shape to the radial ve- 
locity curve in disagreement with the observations. 

In Figure HTH we show the simulated bisector span and 
velocity variations resulting from the blend scenarios de- 
scribed by M0DEL2 and M0DEL3, along with the obser- 
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vations. The expected bisector variations from M0DEL2 
are less than about 25 ms _1 . This is actually consistent 
with the observations, which show no appreciable change 
in the symmetry of the lines within the errors. How- 
ever, the predicted velocity variations from this blend 
scenario are also very small, with a semi-amplitude less 
than 15 ms _1 , which is about 20 times smaller than ob- 
served. MDDEL3 does not fare much better, with a bisector 
semi-amplitude under 50 ms _1 and predicted radial ve- 
locities less than ±30 ms" 1 . Thus, both of these models 
clearly disagree with the observations. 

5. DISCUSSION AND CONCLUDING REMARKS 

The fundamental reason for the failure of all of the 
blend scenarios to reproduce the observations for OGLE- 
TR-56 is that the variations expected in the radial ve- 
locities and in the bisector spans from our simulations 
are of the same order of magnitude, and typically differ 
by a factor of two or less, at least for configurations that 
are consistent with the light curve constr aints. Similar 
results were reported bv lSantos et alJ <)2002|) . In our case 
the predicted bisector span variations tend to be larger 
than the velocity variations. By contrast, the bisector 
variations actually measured for OGLE-TR-56 are essen- 
tially negligible while the velocity variations are highly 
significant. Thus, no model can reproduce both simul- 
taneously, and our tests show this discrepancy cannot 
be solved by adjusting the blend parameters. A further 
reason blend models fail here is the strong observational 
constraint on the brightness of star 1. A star with only 
~3% of the light of star 3 does not produce enough of 
an effect to be noticeable, either in the line profiles or 



in the measured radial velocities. The strong rotational 
broadening of star 1, as predicted from our light curve 
fits, only reduces the effect even further. The conclusion 
must be, then, that a blend configuration involving an 
eclipsing binary along the same line of sight (whether 
in a physical triple system, or in the background) cannot 
account for the observational signatures in OGLE-TR-56 
(photometric and spectroscopic), and we are left with a 
planetary companion in orbit as the only viable explana- 
tion. We are unable to devise a more contrived scenario, 
perhaps involving more stars, that does not go against 
one or another of the observational constraints. 

Dozens of transit candidates are being produced by 
current deep narrow-field searches that focus on very 
faint stars, where crowding is potentially a serious issue 
due to the high stellar densities. Wide-field searches of 
brighter objects have also entered production mode, and 
suffer from a similar confusion problem because of the 
large angular pixel scales. The incidence of blends with 
a background eclipsing binary in both types of surveys is 
therefore expected to be rather high. The same is true in 
most cases for future transit searches from space such as 
NASA's Kepler mission, which will look for even smaller 
photometric signals produced by Earth-size planets. 

The procedures developed in Paper I to model blend 
light curves in detail, along with those described in 21 
to simulate spectral line asymmetries and velocity varia- 
tions, are powerful analysis tools that can be of great 
help for validating these transiting planet candidates. 
The combination of the two techniques makes it possible 
to construct self-consistent and realistic blend scenarios 
and to make predictions that can then be tested quan- 
titatively against the observations. The application to 
OGLE-TR-56 in this paper has provided crucial evidence 
of its true nature, and serves as an example of follow-up 
studies that may be applied to other transit candidates. 
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